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Abstract 

We have carried out deep (V~21 mag) [75 Vi?/ photometric study of the star 
cluster Stock 18. These along with archival Infrared data have been used 
to derive the basic cluster parameters and also to study the star formation 
processes in and around the cluster region. The distance to the cluster is 
derived as 2.8±0.2 kpc while its age is estimated as 6.0 ± 2.0 Myr. Present 
study indicates that interstellar reddening is normal in the direction of the 
cluster. The mass function slope is found to be -1.37±0.27 for the mass 
range 1< M/Mq <11.9. There is no evidence found for the effect of mass 
segregation in main-sequence stars of the cluster. A young stellar population 
with age between 1-2 Myr have been found in and around the cluster region. 
The presence of IRAS and AKARI sources with MSX intensity map also 
show the youth of the Sh2-170 region. 
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1. Introduction 

Stock 18 is a young open cluster and situated in the vicinity of the Sharp- 
less region 170 (Sh2-170), which is almost circular H II region in Cassiopia 

o o 

with angular diameter of about 18 (/ =117.62 , h =+2.27 ; Roger et al. 
2004). The nebula is excited by a centrally located main-sequence (MS) 09- 
type star, BD+63 2093p (Russeil, Adami & Georgelin 2007) which is located 
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at a photometric distance of 2990 pc (Mayer & Macdk 1973). However, deter- 
mination of cluster parameters as well as its photometric study are lacking, 
which are very important to understand the association between the massive 
star, BD+63 2093p, and the cluster Stock 18. In this paper, a comprehensive 
exploration of multi-wavelength data is presented to understand the global 
scenario of star formation in the Stock 18 region. Deep optical (V~21 mag) 
UBVRI photometric data for 13' x 13' region about the cluster Stock 18 
is being reported for the first time in the present study. These along with 
the multiwavelength data are described in the next section while the cluster 
parameters are determined in the section 3. The results obtained from NIR 
data are discussed in section 4 while last section summarizes our results. 

2. Observations and data reduction 
2.1. Optical data 

The CCD fZBVT?/ photometric data were acquired using 104-cm Sam- 
purnanand Telescope (ST) of ARIES, Nainital. The broad-band UBVRI 
observations were standardized by observing stars in the SA 98 field (Lan- 
dolt 1992). In Figure [T]the observed cluster and field regions (see §3.ip are 
marked by circle and rectangle, respectively. The massive stars identified by 
Skiff (2009) along with BD+63 2092p are also marked in Figure [U The log 
of the observations is given in Table [U 

The CCD data frames were reduced using computing facilities available 
at ARIES, Nainital. Initial processing of data frames were done using Iraf 
and Eso-MlDAS data reduction packages. Photometry of cleaned frames 
was carried out using Daophot II software (Stetson 1987). Calibration 
of the instrumental magnitude to the standard system was done by using 
procedures outlined by Stetson (1992). 

For translating the instrumental magnitude to the standard magnitude, 
the calibration equations are derived linear least-square fitting. These equa- 
tions are given below: 

u = U+ (6.941 ± 0.007) - (0.011 ± 0.006)([/ - B) + (0.516 ± 0.012)X 
h = B+ (4.749 ± 0.007) - (0.044 ± 0.005)(5 -V) + (0.274 ± 0.009)X 
v = V + (4.313 ± 0.007) - (0.039 ± 0.005)(F - /) + (0.157 ± 0.009)X 
r = R+ (4.222 ± 0.006) - (0.054 ± 0.007)(V - R) + (0.104 ± 0.007)X 
i = I + (4.711 ± 0.006) - (0.057 ± 0.004)(\/ - /) + (0.047 ± 0.007)X 
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where [/, B, V, R and / are the standard magnitudes, u, b, v, r and i are the 
instrumental aperture magnitudes normahzed for 1 second of exposure time, 
and X is the airmass. The standard deviations of the differences between the 
magnitudes cahbrated using above transformation equations and the stan- 
dards given by Landolt (1992) are 0.035, 0.026, 0.025, 0.026 and 0.027 for 
U, B, V, R and I magnitudes, respectively. The typical Daophot errors in 
magnitude are found to be large (>0.1 mag) for stars fainter than V ~21 
mag, so the measurements beyond this magnitude are not considered. At V 
band, we could detect 2261 stars in nearly 13' x 13' observed region and their 
photometric magnitudes are given in Table 

For further study, it is very important to take into account the incom- 
pleteness that may occur for various reasons (e.g. crowding of the stars). We 
used the Addstar routine of Daophot II to determine the completeness 
factor (CF). The procedures have been outlined in detail in earlier works 
(Sagar & Richtler 1991, Pandey et al. 2001, 2005). The CF values are ob- 
tained to be 100% for brighter stars (V<16 mag). As expected, it decreases 
with faintness of the star and found to be nearly 80% at V~ 21 mag for both 
cluster region as well as field regions (see §3.ip . 

2.2. Archival data sets 

We have also used near- infrared (NIR), mid- infrared (MIR) and far- 
infrared (FIR) archival datasets from the surveys such as Two micron all sky 
survey (2MASS), Infrared astronomical satellite (IRAS), Midcourse space ex- 
periment (MSX) and a Japanese infrared space mission AKARI in the present 
study. FIR data is essential to obtain an unbiased view of star formation and 
in particular to establish whether very young protostars are present (Zavagno 
et al. 2010), especially those undetected at shorter wavelengths. They rep- 
resent the youth of the star forming regions and provide insight into the his- 
tory of star formation. Three IRAS point sources namely IRAS 00001+6417, 
IRAS 23586+6412 and IRAS 23589+6421, and ten AKARI sources namely 
0001312+643829 , 0001577+643734 , 0001275+644107 , 0001039+643907 , 
0001133+643135 , 0000424+644218 , 0002409+643405 , 0000523+643026 , 
0001174+644607 and 0002530+644309 have been identified in the direction 
of cluster. These sources have been detected at least in three bands. AKARI 
and IRAS sources are different, however IRAS source IRASOOOOl+6417 is 
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found within search radius 3" of the position of AKARI source 0002409+643405 
and therefore, they may represent the same source. The presence of these 
infrared sources shows that the surroundings of massive star BD+63 2093p 
provide sites of ongoing star formation which is useful to study the star for- 
mation processes. 



3. Cluster Parameters 

3.1. Radial density profile 

The cluster center is estimated iteratively by calculating average X and 
Y positions of stars with V < 18.0 mag within 80 pixels from eye estimated 
center, until it converged to a constant value. The coordinates of the cluster 
center are found to be RAjsooo = 00^01°^35'^ and DECjaooo = 64°37'37". Us- 
ing this method, a typical error expected in locating the center is 5". The 
center of the cluster coincides (within a limits) with the position of the mas- 
sive star BD+63 2093p, as given by Mayer & Macdk (1973). Therefore, the 
massive star, BD+63 2093p is found to be exactly at the center of the cluster. 

To determine the radial surface density, p(r), we divided the cluster into 
a number of concentric circles. Projected radial stellar density in each con- 
centric annulus was obtained by dividing the number of stars in each annulus 
by its area and the same are plotted in Figure [2] for various magnitude levels. 
The error bars are derived assuming that the number of stars in a concentric 
annulus follow the Poisson statistics. 

The p(r) is parameterized as given by King (1962) 

1 + [r/rc) 

where and /o are the core radius of the cluster and central star density, 
respectively. The field density, /b, is derived using outer region (r > 5'.6) of 
the cluster and 3a of the field density is shown by dashed line in Figure |2j 
We fit the function given by King (1962) to the observed radial density 
distribution of stars by Levenberg-Marquardt non-linear fitting routine and 
plotted by solid lines in Figure HI The best-fit parameters are also given 
in Figure [2l The core radius, Tc, is obtained for various magnitude levels 
(17<V<20), however errors are large in the magnitude range 18 and 20. The 
cluster radius (rd) is thus estimated to be 3f50±0f25, after considering the 
radius at which cluster density is found to be 3a level above the field star 
density, marked by dotted lines in Figure [2J The estimated angular diameter. 
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i.e. 5.7 pc, of the cluster is therefore nearly three times smaller than that of 
Sharpless region Sh-170 as given by Roger et al. (2004). 



3.2. Interstellar extinction 

The interstellar extinction and the ratio of total-to-selective extinction 
Ry = Ay / E{B — V) towards the cluster are important quantities for accurate 
determination of the photometeric distance of the cluster. The extinction 
towards the cluster region is estimated using the {B — V),{U — B) color- 
color (CC) diagram shown in Figure [3l where the MS (Schmidt-Kaler 1982) 
is shifted along the reddening vector having an adopted slope of E{U — 
B)/E{B — V) = 0.72 to match the observations. It shows a spread in the 
observed sequence E(B - V ) = 0.7 to 0.9 mag indicating the presence of 
differential reddening which is an indication of the youth of the cluster. 

To study the nature of the extinction law in the cluster region, we used 
{y — X) vs. {B — V) diagrams, where A is one of the wavelengths of the broad- 
band filters (R, I, J, H, K). This provides an effective method for separating 
the influence of the normal extinction produced by the diffuse interstellar 
medium from that of the abnormal extinction arising within regions having 
a peculiar distribution of dust sizes (cf. Chini & Wargau 1990, Pandey et 
al. 2000). The slope of the distributions rriciuster for colors (V — I), (V — J), 
(V-H) and (V-K) with respect to (B-V) color are 1.11±0.04, 2.08±0.07, 
2.48±0.09 and 2.50±0.10, respectively. These values are similar to that of 
normal reddening towards the direction of cluster. 

The nature of interstellar reddening towards the cluster direction have 
also been analyzed using the color excess ratio method as described by 
Johnson (1968). The intrinsic colors of the stars with spectral types ear- 
lier than AO (V < 14.4 mag) have been obtained using Q-method (Johnson 
&; Morgan 1953) and iteratively estimated the reddening free parameter Q 
[= {U - B) - X{B - V), where X = E{U - B)/E{B - V)], till the color 
excess ratios become constant within the photometric errors (cf. Joshi et 
al. 2008). The color excesses are determined using color relation given by 
Caldwell et al. (1993) for (U - B), (B -V), {V - R) and {V - I) colors 
and by Koornneef (1983) for (V - J), {V - H) and {V - K) colors. The 

, r 1 E{U-B) E{V-R) E{V-I) E(U-B) E{B-V) 

mean vames oi coior excess ratios -^j^—^, ^^^_y^ , ^(^y_j^ , ^{v-j) ' 

If^' ffcjl' ffc?T ffc^ 0.77±0.06, 0.54±0.03, 1.28±0.09, 
0.37±0.04, 0.49±0.04, 0.26±0.02, 0.62±0.02, 1.21±0.03 and 1.26±0.04, re- 
spectively. All color excess ratios are in agreement within la limit that 
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expected for normal interstellar matter (Cardelli, Clayton & Mathis 1989). 
The ratio of total-to-selective extinction is estimated by using the relation 
Rv = I.IE{V - K)/E{B - V) (Whittet & van Breda 1980) and it is found 
to be 2.9 ± 0.3, which is in agreement with the normal value of 3.1 within la 
limit, thus, indicating a normal grain size in the direction of cluster. 

3.3. Optical colour magnitude diagram (CMD) 

The V, {B — V) CMD for the stars in cluster region is shown in Figure HI 
The brightest star BD-l-63 2093p with spectral type 09 is marked by triangle 
while other two early B-type stars Sh2-170 3 and LS 1+64 9 present in our 
photometry are marked by open circles. The UBV magnitudes of the star 
BD-l-63 2093p are taken from Mayer & Macdk (1973), as it gets saturated 
even on our shortest exposure. Using the minimum value of E[B — V) = 
0.70 mag and following relations E{U - B)/ E{B -V)=0.72 and Ay = 3.1 x 
E[B — V), we visually fit the theoretical isochrones (Girardi et al. 2002), for 
log(age) = 6.6 (4 Myr) and 6.9 (8 Myr) with Z = 0.19 to the blue envelope of 
the observed MS. Best-fit isochrone yields a distance modulus of (m — M)v 
= 14.44±0.15 may corresponding to a distance of 2800±200 pc and age=6±2 
Myr. This distance estimation is in agreement with the distance to the star 
BD+63 2093p (Mayer & Macdk 1973). 

Due to better quantum efficiency of CCD detector in the red region, 
y, (y — I) CMD is about a magnitude deeper in comparison to V, {B — V) 
CMD. Hence, V, (V^ — /) CMD is useful to probe the stars at lower mass range 
and the same has done below. Figures |5] (a) and (b) show V, {V — I) CMDs 
for the cluster region and field region, respectively. The contamination due 
to field population is clearly visible in the CMD of cluster region. To remove 
contamination of field stars from the MS and PMS sample, we statistically 
subtracted the contribution of field stars from the CMD of the cluster region 
using the procedure called zapping described in Sandhu, Pandey & Sagar 
(2003). For a randomly selected star in the V, {V - I) CMD of the field 
region, the nearest star in the cluster's V, {V — I) CMD within V ± 0.35 
and (V^ — /)±0.2 of the field star was removed. While removing stars from the 
cluster CMD, necessary corrections for incompleteness of the data samples 
were taken into account. The statistically cleaned V, {V — I) CMD of the 
cluster region is shown in Figure [5] (c) which clearly shows the presence of 
PMS stars in the cluster. 

Figure |6] shows statistically cleaned unreddened Vq, {V — /)o CMD. The 
stars having spectral type earlier than AO were individually unreddened (cf. 



6 



§3.2p . whereas mean reddening of the region, estimated from available indi- 
vidual reddening values in that region, was used for other stars (cf. Sharma 
et al. 2007 ). The isochrone for 4 Myr by Girardi et al. (2002) and PMS 
isochrones by Siess, Dufour & Forestini (2000) have been plotted in Figure [61 
Evolutionary tracks by Siess, Dufour & Forestini (2000) for various masses 
have also been shown in the same figure. 

3.4- Mass function 

With the help of statistically cleaned CMD, shown in Figure [6l we can 
derive the mass function (MF) using the theoretical evolutionary model of 
Girardi et al. (2002). Here, we are not discussing the mass function of the 
stars brighter than Vq < 9.5 mag (V^ <12.0 mag) because they are saturated 
even on our shortest exposure. Since post-main-sequence age of the cluster 
is 6 Myr, the stars having Vq < 13.5 (M >2.5M0) have been considered to 
be on the MS. For the MS stars, the LF was converted to the MF using the 
theoretical model by Girardi et al. (2002) (cf. Pandey et al. 2001, 2005). 
The MF for PMS stars was obtained by counting the number of stars in 
various mass bins (shown as evolutionary tracks; Siess, Dufour & Forestini 
(2000)) in Figure ini The MF for the whole cluster is given in Table |3] and 
plotted in Figure [71 The slope of the MF for the entire observed mass range 
1.0< M/Mq <11.9 is derived to be F = -1.37 ± 0.27, which is in agreement 
with the Salpeter value ( = -1.35; Salpeter 1955). 

3.5. Mass segregation 

To characterize the degree of mass segregation in Stock 18, the MS sample 
is subdivided into two mass groups (5.55 < M/M© < 11.88, 2.43 < M/M© < 
5.55). Figure [HI shows cumulative distribution of MS stars as a function of 
radius in two different mass groups. The figure indicates that more massive 
stars (5.55 < M/M© < 11.88) tend to lie toward the cluster center. How- 
ever, the Kolmogorov-Smirnov test shows that the confidence level with the 
statement is very low, i.e., 45%. Therefore, no statistically significant effect 
of mass segregation is found in the cumulative distribution. 

4. 2MASS NIR data 

The dust grains in the circumstellar disk of PMS stars absorb a fraction 
of visible or ultraviolet photons emitted by central star, and re-radiates the 
absorbed energy in IR wavelengths. Thus, the presence of NIR excess can be 
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used to distinguish between the PMS stars from other stars in optical CMD. 
Therefore, NIR observations are very effective for investigating the nature 
of obscured clusters and the populations of YSOs which are embedded in 
molecular clouds (Lada & Adams 1992). 

Figure [9] shows (J — H), {H — K) CC diagram for the cluster region 
(r < Tci) and 0.5 degree region centered around the cluster which represents 
the Sharp less region Sh-2 170. Unreddened MS and giant branch, taken from 
Bessell & Brett(1988), are represented by solid lines. The reddening vectors 
for early and late-type stars are shown by parallel solid lines drawn from the 
base and the tip of two branches. Location of T-Tauri stars (marked by TTS; 
Meyer, Calvet & Hillenbrand 1997), proto-star (marked by PS) like objects 
and Herbig Ae/Be (Hernandez et al. 2005) are also shown. The extinction 
ratio Aj/Av = 0.282, Ah/Av = 0.180 and Ak/Av = 0.116 have been taken 
from Cardelli, Clayton & Mathis (1989). Stars below the reddening vectors 
in TTS region in Figure |9] are considered to be NIR excess stars. These 
stars are young with the characteristics of having circumstellar material. 
We found six such stars in the cluster region, however, the Sharpless region 
Sh-2 170 contains eleven NIR excess stars. Therefore, nearly 54% of total 
stars present in Sharpless region Sh-2 170 are found within cluster region 
(r < Tci). Carpenter (2000) estimated that 55% of the YSOs present in SFRs 
are contained in the clusters. However, Bressert et al. (2010) showed that 
only a low fraction (<26%) of YSOs are formed in dense environments using 
Spitzer Space Telescope surveys. Out of the six NIR excess stars in cluster 
region, only four stars have optical counterpart. 

In Figure [TOl spatial distribution of 0/B-type stars (plus symbol), IR- 
excess sources (open circles), IRAS points sources (crosses) and AKARI 

sources (open squares) are displayed on a 30 x 30 DSS-II R band im- 
age around the cluster. The MSX A-band intensity map is superimposed 
on DSS-II R band image in Figure [TOl which indicates presence of several 
discrete sources representing high-density clumps and region of further star 
formation. The ^^CO temperature map of the region has also been shown in 
the Figure [HJ] using the data from Kerton & Brunt (2003). The AKARI and 
IRAS sources are associated with MSX intensity map as well as with ^^CO 
map, indicate the youth of these regions. The presence of FIR IRAS and 
AKARI sources which may even younger than 1 Myr, supports the ongoing 
star formation process within Sh-2 170 region. 

The locations of NIR excess sources within cluster region in optical Vq, {V— 
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/)o CMD are shown in Figure [TTl The isochrone for 4 Myr by Girardi et al. 
(2002) and PMS isochrones by Siess, Dufour & Forestini (2000) are also over 
plotted in Figure [TTJ These probable YSOs are individually unreddened us- 
ing NIR data. The NIR CC diagrams were used to estimate Ay for each of 
these YSOs by tracing them back to the intrinsic CTTS locus of Meyer et 
al. (1997) along the reddening vector (for details see Ogura et al. 2007). It 
appears that the age of these NIR excess PMS stars is 1-2 Myr. However, the 
estimation of the age of the PMS stars by comparing the observations with 
the theoretical isochrones is prone to the random as well as systematic errors 
(see e.g., Hillenbrand 2005, Hillenbrand et al. 2008). The other sources of 
errors which can affect the estimation of ages may be the intrinsic variability 
of YSOs, accretion processes, uncertainties in determination of extinction, 
binarity and dispersion in distance measurement (Hartmann 2001). Binarity 
will brighten the star, consequently the CMD will yield a lower age estimate. 
In the case of equal mass binaries, we expect an error of 50 - 60% in the age 
estimation of the PMS stars. However, it is difficult to estimate the influ- 
ence of binarity on the mean age estimation as the fraction of binaries is not 
known. If the PMS stars are cluster members, it appears that there is differ- 
ence in age of YSOs and massive star BD +63 2093p. Moreover, the PMS 
stars are located in the clumps (please see the Figure 10) and star BD+63 
2093 is located out side of these clumps within the cluster region. In past, 
it was also noticed that ongoing star formation may occur in the vicinity of 
massive star with in a very small angular diameter even 1 pc ( e.g. Azimlu 
& Fich 2011). However, current data are not enough to make a confident 
claim on the ongoing star formation and more FIR and deep optical data is 
needed to probe the real picture of star formation in the region of the cluster 
Stock 18. 

5. Summary and Conclusion 

On the basis of a comprehensive multi- wavelength study, here, an attempt 
has been made to understand the basic parameters and the implications of 
star formation processes in the cluster Stock 18. Deep optical UBVRI data. 
along with archival data from the surveys such as 2MASS, MSX, IRAS, 
AKARI and ^^CO are used to understand the global scenario of star forma- 
tion in and around the cluster region. Reddening, E{B — V), in the direction 
of cluster is found to be normal but varying between 0.7 to 0.9 mag. The 
post- main- sequence age and distance to the cluster are found to be 6±2 Myr 
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and 2.8±0.2 kpc respectively. The slope of the MF is found to be -1.37±0.27 
for the mass range 1.0< M/Mq <11.9. The effect of mass segregation is 
not seen in the MS stars. A population of PMS YSOs in the cluster region 
having masses 0.8 - 2 Mq have been found within the cluster. The position 
of the YSOs on the CMDs indicates that these stars have age between 1 to 2 
Myr. However, this age spread is affected with observational errors and lim- 
itations of theoretical models. Further, the present data sets are not enough 
to explain the ongoing star formation. 
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Table 1: Observation log of CCD observations of the cluster Stock 18 and the calibration 
region SA98 (Landolt 1992). 



Date(UT) Filter Exposure Time (s) 

(xno. of exposures) 







Stock 18 


SA 98 


08 October 2007 


U 


1800x2 






B 


1200x3 






V 


900x3 






R 


300x3 






I 


300x3 




19 December 2006 


V 


900x3 






R 


480x2,60x1 






I 


480x2,300x2 




23 December 2006 


U 


300x2 


300x7 




B 


120x3 


120x7 




V 


120x3 


60x7 




R 


40x3 


40x6 




I 


50x5 


10x7 



Table 2: UBVRI photometric data of the sample stars in 13'xl3' region. 
CCD positions of stars are converted into RAj2ooo and DECj2ooo using the 
Guide Star Catalogue II (CSC 2.2, 2001). The complete table is available in 



electronic form. 



ID 


R'Aj2000 


DEC,J2000 


U 


B 


V 


R 


/ 




(deg) 


(deg) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


0.388444 


64.612169 


11.639 ± 0.008 


12.016 ± 0.002 


11.608 ± 0.008 


11.386 ± 0.002 


11.089 ± 0.002 


2 


0.381583 


64.677193 


12.609 ± 0.008 


12.878 ± 0.003 


12.241 ± 0.004 


11.882 ± 0.007 


11.453 ± 0.006 


3 


0.417500 


64.531441 


14.559 ± 0.038 


13.639 ± 0.005 


12.247 ± 0.004 


11.436 ± 0.026 


10.595 ± 0.031 


4 


0.306694 


64.588333 


12.618 ± 0.009 


13.007 ± 0.004 


12.542 ± 0.003 


12.274 ± 0.003 


11.938 ± 0.003 
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Table 3: The MF of the cluster Stock 18. The numbers of probable cluster members (N) 
have been obtained after subtracting the expected contribution of field stars in i)3.3l log^ 
represents log {N/d{log m)). 



Magnitude range 


Mass Range 


Mean 


Whole region 


(Vo mag) 


(Mo) 


log (Mq) 


N 


log0 




9.5 - 10.5 


11.88 - 8.17 


1.001 


4 


1.391 


10.5 - 11.5 


8.17- 5.55 


0.836 


2 


1.076 


11.5 - 12.5 


5.55 - 3.52 


0.657 


8 


1.607 


12.5 - 13.5 


3.52 - 2.43 


0.473 


9 


1.748 




2.50 - 1.50 


0.301 


27 


2.085 




1.50 - 1.00 


0.097 


52 


2.470 
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0:02:00.0 30.0 
RA(J2a00) 



01 :00.0 



Figure 1: Observed image of Stock 18 from ST Nainital. The observed cluster and field 
regions are marked as circle and rectangle, respectively. Position of the star BD+63 2093p 
and other massive stars from Skiff (2009) are marked by the symbols of cross and plus, 
respectively. 
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Figure 2: Projected radial stellar density profile of Stock 18. Dashed lines represent 3cr 
levels above the field star density and solid curve shows best fit to the empirical model of 
King (1962) . 
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Figure 3: The (B — V), {U — B) CC diagram for stars in the chister region. Sohd hnes are 
intrinsic MS reddened along the reddening hne with E{B — V) ~ 0.7 mag and 0.9 mag. 
BD+63 2093p is marked by triangle, and early B-type stars Sh 2-170 3 and LS I +64 9 
present in our photometry are marked by open circles. 



17 




(B-V) 

Figure 4: V, {B — V) CMD for the cluster regfon. Theoretical isochrones from Girardi et 
al. (2002) are shown with continuous solid and dotted lines for log(age)= 6.6 and 6.9 Myr, 
respectively. Early B-type stars Sh 2-170 3 and LS I +64 9 are present in our photometry 
and marked by open circles. Massive star BD+63 2093p is marked by triangle and its 
UBV data is taken from Mayer & Macdk (1973). 
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Figure 5: V, {V — 1) CMDs for (a) cluster region (b) field region (c) the statistically 
cleaned CMD for cluster. Theoretical isochrones from Girardi ct al. (2002) are shown 
with continuous and dotted lines for log(age)= 6.6 and 6.9 Myr, respectively. Early B- 
type stars Sh 2-170 3 and LS I +64 9 are marked by open circles. 
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Figure 6: Statistically cleaned Vo/{V — I)q CMD for stars lying in the cluster region. The 
isochrone for 4 Myr age by Girardi et al. (2002) and PMS isochrones of 1,2,5,10 Myr along 
with evolutionary tracks of different mass stars by Siess, Dufour & Forestini (2000) are 
also shown. All the isochrones are corrected for a distance of 2.8 kpc. Early B-type stars 
Sh 2-170 3 and LS I +64 9 are marked by the symbol of asterisk. 
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Figure 7: A plot of the MF in the cluster, log(^) represents log {N/d{log m)). The error 
bars represent ±\/]V errors. The solid line shows a least square fit to the entire mass range 
1.00 < M/Mq < 11.88. 



21 




22 




0.6 1 1.5 
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Figure 9: CC diagrams using the 2MASS JHK data, upper panel: The NIR CC diagram 
for the cluster region Stock 18. lower panel : Same as the upper panel, but for the 0.5 
degree surrounding region centered around the cluster and represent the Sharpless region 
Sh2-170. Filled circles in TTS region denote T-Tauri stars with NIR excess. 
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03:00,0 02:00.0 01:00,0 0:00:00.0 
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Figure 10: Spatial distributions of 0/B-type stars (plus symbol), IR-excess sources (prob- 
able T-Tauri type stars, small open circles), IRAS points sources (crosses) and AKARI 
sources (open squares) are overlaid on the image of Sh2-170 in DSS-II R band. Cluster 
region Stock 18 (r < rd) is marked by big open circle. ^^CO contours (black contours) 
from Kerton & Brunt (2003) and MSX A-Band intensity contours (white contours ) have 
also been shown. The MSX A-Band contours are at 0.12, 2.24, 2.41, 2.55, 2.68, 2.85, 5.26 
xlO-'^ W m-2 Sr-i. ^^CO contours are at 9.0, 9.3, 9.6, 9.9, 10.2, 10.5 K. 
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Figure 11: NIR excess sources (see 21) within the cluster region in Vq/{V — I)o CMD. 
The isochrone for 4 Myr age by Girardi et al. (2002) and PMS isochrones of 1,2,5,10 Myr 
along with evolutionary tracks of different mass stars by Siess, Dufour & Forestini (2000) 
are also shown. All the isochrones are corrected for a distance of 2.8 kpc. 
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